Pedestrian Flow Models with Slowdown Interactions 

Alina Chertock* Alexander Kurganov} Anthony Polizzif and Ilya Timofeyev§ 



Abstract 

In this paper, we introduce and study one-dimensional models for the behavior of pedes- 
trians in a narrow street or corridor. We begin at the microscopic level by formulating a 
stochastic cellular automata model with explicit rules for pedestrians moving in two oppo- 
site directions. Coarse-grained mesoscopic and macroscopic analogs are derived leading to 
the coupled system of PDEs for the density of the pedestrian traffic. The obtained PDE 
system is of a mixed hyperbolic-elliptic type and therefore, we rigorously derive higher-order 
nonlinear diffusive corrections for the macroscopic PDE model. We perform numerical ex- 
periments, which compare and contrast the behavior of the microscopic stochastic model and 
the resulting coarse-grained PDEs for various parameter settings and initial conditions. We 
also demonstrate that the nonlinear diffusion is essential for reproducing the behavior of the 
stochastic system in the nonhyperbolic regime. 

1 Introduction 

In contrast with the considerable effort devoted to the modeling of vehicular traffic, modeling of 
pedestrian traffic received little attention until fairly recently. In the past two decades, a variety 
of pedestrian traffic and crowd dynamics models have been proposed, ranging from agent-based 
microscopic to macroscopic PDE models describing various phenomena such as crowd behavior 
under panic, pedestrian planning, structural design, etc., see [3|IT2|[T4"] and the review papers [^[13] 
for examples and references. 

Cellular automata (CA) has a long history of applications in different areas of science and 
engineering. In particular, the CA models have been applied to the vehicular traffic to derive a 
coarse-grained PDE description with the look-ahead dynamics for the density of the car traffic [20] . 
Multilane and multiclass were also considered as extensions of the original model [TJ[7JEJ|21] • CA 
microscopic models with empirical rules have also been used to simulate pedestrian movement 

psisiin]. 

The main emphasis of the present work is on the connection between the microscopic CA models 
for the bi-directional pedestrian traffic and its coarse-grained PDE analogs. In particular, we 
utilize the CA approach to formulate a one-dimensional (1-D) microscopic model for the pedestrian 
motion in a narrow street of corridor and derive a corresponding PDE description for the density 
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of the pedestrian flow. The present work is motivated by the recent pedestrian experiments and 
modeling discussed in [2j[T8], where the authors consider the pedestrian motion in a circle and 
model the crowd movement using a system of 1-D PDEs. 

The major advantage of the CA formalism is that it allows for a systematic derivation of the 
coarse-grained dynamics. The main assumptions about the traffic flow (vehicular or pedestrian) 
are build into the microscopic model. The derivation of the coarse-grained description typically 
requires some simplifying assumptions about the statistical behavior of the microscopic model 
which can be numerically verified. 

The main conceptual difference between the CA models of the vehicular traffic and the CA 
model studied in this paper is that the pedestrian motion is bi-directional which leads to completely 
new paradigms. In particular, the coarse-grained description of the pedestrian traffic becomes a 
system of conservation laws of a mixed hyperbolic-elliptic type. We demonstrate that this system 
may exhibit an unrealistic nonhyperbolic behavior, depending on the magnitude of the density 
of the pedestrian traffic. To overcome this difficulty, we systematically derive nonlinear diffusive 
corrections to the PDE model from the original microscopic description. 

The paper is organized as follows. In Section |2j we introduce the microscopic CA model. In 
Sections [3] and HJ we discuss the derivation of the mesoscopic and macroscopic analogs for the 
density of the pedestrian traffic, as well as derive the next-order nonlinear diffusive corrections 
for the macroscopic PDE model. In Section [5], we describe several sets of simulations, in which 
the microscopic and macroscopic PDE models are compared illustrating, in particular, the effect 
of diffusion for the initial conditions leading to the regime of nonhyperbolicity of the inviscid 
macroscopic PDE model. Finally, in Appendix A, we provide a brief description of the numerical 
method used to solve the derived systems of PDEs. 

2 Microscopic Pedestrian Model 

To construct the microscopic model, we consider the time evolution of a 1-D lattice L with pedes- 
trians moving into two opposite directions. The process is schematically illustrated in Figure 12.11 
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Figure 2.1: Schematic pedestrian configuration: Filled and empty circles represent pedestrians moving 
to the right and left, respectively. Arrows represent possible pedestrian movements in this configuration. 



The main difference with previous works on car traffic models (e.g., [20]) is that the pedestrians 
can move into opposite directions and we need to specify explicit rules when two such pedestrians 
meet. We represent pedestrians moving to the right and to the left by the variables &^{t) = cr + (k, t) 
and cr^(t) = a~(k,t), respectively, where k is the position on the lattice and t denotes time: 



1, if at time t cell k is occupied by a pedestrian moving to the right (left), 
0, otherwise. 
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Moreover, we assume that pedestrians moving into two opposite directions can occupy the same 
cell (i.e., it is possible to have (J k (t) = cr k (t) = 1 for some k), but their velocities are reduced 
in that case. This mimics the slowdown due to side-stepping in more realistic CA models [5l[9] 
when pedestrians are avoiding the upcoming traffic. We also assume that pedestrians moving in 
the same direction cannot occupy the same cell. For simplicity of the presentation, we omit the 
notation for the time-dependence of a k and a k for the rest of the paper. 

To complete the construction of the microscopic CA model, we consider explicit rules for the 
slowdown interaction which describe how the velocity changes when pedestrians moving in opposite 
directions come in contact (i.e., occupy the same cell or two neighboring cells). In particular, we 
explicitly prescribe velocities for four different pedestrian configurations in the cells neighboring 
to the right-moving pedestrian with a k = 1 (assuming that cr k+1 = 0): 

' Co, if a k = cr k+1 = (no left-moving pedestrians in cells k or k + 1), 

ci, if <7fc = 1, cr A T +1 = (a left-moving pedestrian is in cell k), 

c 2 , if a k = 0, cr A T +1 = 1 (a left-moving pedestrian is in cell k + 1), 

k C3, if a k = cr k+1 = 1 (left-moving pedestrians in cells k and k + 1). 

Notice that from the common sense considerations, the velocities should obey the following rela- 
tionship: C3 < c 2 ~ c\ < Co. The velocities of the left-moving pedestrian a k = 1 can be obtain in 
a similar manner. 

The prescribed velocities are used to determine the probability of a pedestrian to move to the 
neighboring cell, that is, the probability of transition k — > k + 1 for a k during a small time interval 
At is 



P^k+i = ^ [coa+(l - <r+ +1 )(l - - a- +l ) + Cl cr+(1 - a+ +1 )a^l - a k+l ) 

+C 2 (7+(1 - (7+ +1 )(l - <7fc )<7fc +1 + c 3 a+(l - at+i)°k°k+i 
while the probability of transition k — > k — 1 for o k is 

CQCTfc (1 - <7fc_i)(l - 0fc_i)(l - fffc ) + Cl<7fc (1 - <7*_i)(l 

+c 2 a k (l - a^a^l - a+) + c 3 a k (l - v^W-i^k 



(2.1) 



(2.2) 



The lattice configurations er t ~ := {cr k } and erf := {& k } together constitute a continuous-time 
Markov chain for cr t := {cr^ , erf}. This model is easily simulated numerically using the Metropolis 
algorithm for computing the expected values of a~ and a + . One can develop a Kinetic Monte- 
Carlo algorithm for these problems, but we found that the Metropolis algorithm was quite efficient 
when the velocities cq, c\, c 2 , C3 are not very small. 

Since er t is a continuous-time stochastic process, its generator is defined by 

E^(o-At) - 1p(<T ) 

Aw = lim , 

r At^O At 

where <To is the initial configuration, er^t is the configuration at time At, i/j is any test function, 
and the expectation is taken over all possible transitions from er to cr^ t . The generator of this 
stochastic process can be computed as follows: 

^ = Xt £ { P ^+i W ff iU+i> °o) - lK<*>)] + P^k-x , - ^o)] } , (2.3) 
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where cr^ fc+1 is the configuration obtained from cr^ by exchanging the values in cells k and k + 1, 
and similarly, cr^ k _ 1 is the configuration obtained from cr^ by exchanging the values in cells 
k and k — 1. To derive the coarse-grain model, we need to compute the generator (12. 3p in two 
particular cases: For ^(cr + , <x~) = and V ; ( cr+ > °"~) = &k > we obtain 

, + = ^^w and A n+i-^fc-^fc-i (24) 

fc At fe At v y 

respectively. 



3 Mesoscopic Model 

In this section, we use the microscopic CA model presented in Section [2] to derive the mesoscopic 
model for the densities p^ := Ecr^ (once again, to simplify notation, the time-dependence of p^ 
is omitted throughout the paper). To this end, we first recall that the generator A satisfies the 
following property: 

= EAip, 

at 

which can be applied to the test functions -0 = a£ and tp = This together with ( 12. li . ( 12. 2ft 
and ( 12. 4p results in the equations for the time-evolution of pf: 



dpi 
dt 



+Clf7+_ 1 (1 - C7+)c7fc_ 1 (l - C7fc ) - Ci(7+(1 - t7+ JcT^ (1 



d Pk 
dt 



(3.1) 



k+l) u k V x ^fc+lJ 
+C 2 ff+„ 1 (1 - <7+)(l - <7fc_!)<7fc - C 2 CT+(l - (T+ +1 )(l - (7^)(7^ +1 

+C 3 0-fc_i(l - ^fc )^fe-lO-fc - C 3 (7+(l - 0"+ +1 )c7fe (7fc + J , 

E[co^ +1 (l - a* )(1 - " <i) " cb<Tfc(l - - - 

+c 2 a fc+1 (l - fffc )(7+(l - (7+ +1 ) - c 2 a k {l - Vk-iW-ii 1 - a t) 



The system ( 13. ip is exact, but not closed, since its right-hand side involves higher-order mo- 
ments. The closure approximation can be derived by assuming that the joint measure for a + 
and a~ is approximately a product measure. This implies, for instance, that adjacent cells are 
approximately independent and, in particular, E[<7^<7^ +1 ] ~ E<7 fe E<7 fe+1 . Thus, a closed system of 
equations for pf can be obtained and the resulting mesoscopic model for the pedestrian density 
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reads 



do + 

-jt =c,pU{\ - P +)(i - p-^Ki - p k ) - c oP i(i - - - Pk+1 ) 

+Cipt-x(l ~ Pk)Pk-li l ~ Pk) - ClPfc(l - Pk+i)Pl (1 - Pfc+l) 
+C2pt-!^ - Pt )(! - Pk-i)Pk ~ (! - Pfc+iX 1 - Pfc K+l 

+C3Pfc-l( 1 - Pk)Pk-lPk - °3pt (1 - Pk+l)PkPk+H 



d Pk 



(3.2) 



=c p fc+1 (l - P fc )(1 - - - c p fc (1 - p fc „ x )(l - p^_ 1 )(l - Pt 



dt 

+ClPfc +1 (l - Pfc )(1 - Pfc )Pfc+i - CiPfc (1 - Pfc_J(l - Pfe_i)Pfc 

+c 2 p^ +1 (i - Pfc )p^(i - p£t-i) - c 2 p^ (i - Pfe_ 1 )p^_ 1 (i - pt) 
+ c 3Pfe+i( 1 - pk)pipi+i - c aPfc (i - Pk-i)pt-iPk ■ 

Note that the system (13. 2p is defined on the same lattice £ as the microscopic model. 

4 Macroscopic PDE Model 

We now treat sites k e L as cells with some fixed length /i > 0. Let Q denote a sub domain of R 
corresponding to the lattice £, i.e., f2 = [0, L] (the number of cells, of course, depends on h). We 
consider a rescaling of time t — > ht and derive a coarse-grained PDE model in the limit as the cell 
size tends to zero and the number of cells tends to infinity. 

To this end, we rewrite the system (13. 2p in the following flux form (taking the time rescaling 
into account): 

dpj _ Fk,k+i ~ -^fc-i,fc dp k _ Fk,k+i ~ Fk-i,k .q 
dt h dt h 

where 

Kk+i = Pki 1 - Pk+i) [C 1 - Pfc+i) M 1 - Pk) + c iPfe) + Pfc+i ( c 2(l - Pfe) + csPfc )] , ^ 2 ^ 

^Tfc+i = Pk+ii 1 - Pk) [(! - Pfc ) ( c o(l - Pfc+i) + CiPi+i) + Pfc (c 2 (l - Pfc+i) + c 3 p^ +1 )] • 

Multiplying the above equations by tpk ■= <f(kh), where ip G Cg(f2) is a test function, and using 
the summation by parts property over Q, yields 

k k 

Next, we multiply equation (14. 3 j) by /i and expand <pfc+i in a Taylor series about fc/i to obtain 

Z>*^f* = ±E F M+M + (4.4) 

We define pedestrian densities on Q as follows. Again, using the notation p* 1 (for convenience), 
define the function p ± (x,t) as a continuous piecewise linear interpolation (in the spatial variable) 



of pi(t) and take the limit as h — > + . Due to the boundeness of both p^ 1 and we obtain a 
weak formulation of the coarse-grained model: 

J <p(x)^- t p ± (x J t)dx = ± j F ± (p + ,p-) V f (x)dx, (4.5) 
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where F ± (p + ,p ) are defined as the corresponding limits of Fj^ k+V i.e., 

F + (p + ,p~) = f(p + )g(p~), F-(p + ,p~) = f(p-)g( P + ), 

where 

f(u)=u(l-u), g(u) = (c 3 - c 2 - ci + c )u 2 + (c 2 + c x - 2c )w + c . (4.6) 
Since <£> is arbitrary, the integral equations H4.5[) can be written as the following system of PDEs: 

B 4 + l IMP-)] = <>, % - °- [«,->,(,♦>] = 0. (4.7) 



Remark 4.1 Note that the velocities ci and c 2 enter only as a sum into (14. 6p . Therefore, it is 
not necessary to specify them separately. 

Remark 4.2 The coarse-grained system ( 14 .71) is only conditionally hyperbolic. Indeed, the Jaco- 
bian 

f'(p + )g(p-) f(p + )g'(p') 
-f(p~)g'(p + ) -f'(p~)g(p + ) 

has real eigenvalues only if 

[f'(p~)g(p + ) + f'(p + )g(p-)} 2 - ±f(p-)f(P + )g\p-)g\p + ) > o. (4.9) 

Therefore, for any particular choice of velocities Cq, c±, c 2 , and C3 there is a region on nonhyperbolic- 
ity in the (p~, p + ) plane. From (14.61) . one can see that the nonhyperbolicty can only manifest itself 
when pedestrians moving in two opposite directions are both present in a particular location. Two 
examples of nonhyperbolic regions are plotted in Figure 14.11 The nonhyperbolic region described 
in (14.91) depends only on the ratio of velocities c\/co, c 2 /co, and cs/c , but not on the particular 
value of Co- The region of nonhyperbolicity becomes larger as the slowdown effect becomes more 
pronounces (i.e., as the ratios ci/c , c 2 /c , and c 3 /c become smaller). The loss of hyperbolicity 
may induce instabilities (as illustrated in Section [5]), which are nonphysical and can be removed by 
introducing a nonlinear diffusive correction to the system (as we demonstrate in the next section). 




4.1 Diffusive Correction 

The derivation of the coarse-grained PDE system 04 .7p . (14 .6p can also be obtained by formally 
using the Taylor expansions 

pLi = pt±Hpty+Y^y +0 ^ 

in the flux formulation (14. ip . (I4.2p followed by passing to the limit as h — > + . 
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Figure 4.1: Regions of nonhyperbolicity (inside) in the p + ,p~ plane for the equation (14. 7p for two 
particular choices of velocities; left part c = 1, c\ — c 2 = c /2, c 3 = c /4 and right part c = 1, 

ci = c 2 = cq/3, c 3 = co/6. 



Alternatively, keeping /i fixed and neglecting the (9(/z 3 ) terms, leads to the following second- 
order PDE system, which contains nonlinear diffusion terms: 



<9p H 







dt dx 



[f(P + )9(P~)] 



c 



ptx + (ci - c + (c 3 - c 2 - Cl + c )p + (c 2 - ci)p + ) p x p+ + 



-(ci - c 2 )p + (l - p + )p x:c + - ((ci + c 2 - 2c )p + (c 3 - c 2 - d + c )(p ) 2 ) pi 



h 



yPxx + (ci - c + (c 3 



c 3 - ci + c )p + + (c 2 - ci)p ) plp x + 



-(ci - c 2 )p (1 - p )p+. + - ((c 2 + ci - 2c )p + + (c 3 - c 2 - ci + c )(p + ) 2 ) p~ 

The diffusion in the above system can be simplified considerably for certain conditions on the 
velocities. In particular, assuming that 

ci = c 2 , (4.10) 

and replacing h (a fixed cell size in the mesoscopic model) with a small parameter e, the diffusive 
correction to the original coarse-grained PDE model (14. 7p becomes 



(4.11) 



¥ + Tx WW-)] = f [((1 - P~f + 2a ip -(l - p~) + a 3 (p-) 2 ) pt] x , 
where «i = Ci/c and a 3 = c 3 /c . 

Remark 4.3 Notice that the assumption (14.101) is rather mild since the velocities c\ and c 2 only 
enter as a sum into the fluxes (see Remark 14.11) . We would also like to stress that the coefficients 
of the nonlinear diffusion are positive provided both p + and p~ are between and 1. 

Remark 4.4 The nonlinear diffusion in (14.111) reflects the presence of pedestrians moving in the 
opposite direction. For instance, the diffusion coefficient in the first, p + -equation depends on p~ . 
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It is instructive to consider the limiting cases: If p~ — (i.e., no pedestrians moving to the left 
are present), the diffusion reduces to the usual linear diffusion 0.5ec p xx , while if p~ = 1, then 
the diffusion becomes 0.5ec o a 3 p xx = 0.5ec 3 p xx . The latter is, again, a linear diffusion, but with a 
smaller coefficient (since c 3 < c ) reflecting a high density presence of the pedestrians moving in 
the opposite direction. 

Remark 4.5 It should be observed that the size of the viscosity coefficient e needs to be estab- 
lished experimentally for each particular problem at hand. In our numerical examples reported in 
Sections 15.21 and 15. 3^ the value of e is chosen empirically. 

5 Simulations 

In this section, we present several sets of numerical experiments comparing and contrasting the 
propagation of pedestrian density in ensemble simulations of the microscopic stochastic model 
and the corresponding numerical solutions of the macroscopic PDE models (with and without 
diffusion). The mesoscopic model ( 13. 2 p can be considered as a conservative first-order finite- 
difference discretization of the PDE models and the results of the mesoscopic simulations (not 
shown in the paper), performed on a sufficiently fine grid, typically agree very well with the PDE 
results. 

In all of the PDE simulations below, we implement a semi-discrete second-order central-upwind 
scheme from [15l[16] with the minmod parameter 8 = 1, CFL number equal to 0.5. The scheme is 
briefly described in Appendix [A] It should be observed that the PDE simulations are performed 
on a much coarser mesh than the corresponding microscopic ones. We use periodic boundary 
conditions in all of the presented numerical examples. 

In the first set of numerical experiments (Section 15. ip . the initial conditions resemble the "red 
light" situation when a group of tightly packed pedestrians with density 1 is released at time 
t = 0. In the second set (Section I5.2p . we mimic the pedestrian movements starting with fully 
mixed initial conditions sampled from a particular piecewise constant density. Finally, in Section 
I5.3[ we illustrate the stabilizing effect of the diffusive corrections in the nonhyperbolic regime. 

5.1 "Red Light" Initial Conditions 

We first consider the microscopic CA model with the following initial conditions: 

a^,0) = { 1 ' ^ fc ^> a-(fc,0) = ( 1 ' N-n 2 <k<N- ni , 
10, otherwise, 10, otherwise, 

with ri2 — rii <C N. These initial conditions correspond to two (relatively small) groups of pedes- 
trians standing still and starting to move toward each other at time t = 0. The velocities are 
taken as 

c = 0.8m/ s, ci = c 2 = c /a, c 3 = c /(2a), (5.2) 

where the parameter a describes the strength of the slowdown interactions, for which we consider 
two regimes with either a = 2 or a = 3. The parameters in the CA model are chosen so that 
the two groups are away from the boundary and, thus, boundary conditions do not affect the 
interaction: 



N = 1400, m = 301, n 2 = 340, h = 0.2m, At = 0.01s, MC = 5000, (5.3) 
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where MC is the number of Monte-Carlo simulations. 

The macroscopic simulations of the purely convective PDE system (14. 7p are performed on the 
computational domain [0, 280] with the mesh size Ax = 0.8 and initial data corresponding to (15. ip : 



p + (x,0) 



60 < x < 68, 
otherwise, 



p-{x,0) 



212 < x < 220, 
otherwise, 



(5.4) 



with the velocities c , Ci, c 2 and c 3 same as in (15. 2p . 

The initial settings (I5.ip - ()5.3p and ( 15 .4|) correspond to the right- and left-moving groups, which 
initially do not overlap and their dynamics is equivalent to moving cars. After some time, the 
two groups begin to interact and we study how well the macroscopic PDE model reproduces 
these interactions. Figures 15.11 and 15.21 show comparison of the density profiles in CA and PDE 
simulations. 



Time = 
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Figure 5.1: Comparison of the density profiles computed by the CA (solid line) and PDE (dash-dotted 
line) models with a = 2 and the "red light" initial conditions in ( 15.11) and ( 15.41) . The right- (p + ) and 
left-moving (p~) pedestrians are represented by the blue and red colors, respectively. 

The macroscopic PDE model ( 14 .7p produces very good results in the a = 2 regime. In par- 
ticular, it reproduces the forward front at times t = 80, 110 and 140 quite well. The macroscopic 
results are slightly different from the stochastic ones after the two groups pass through each other 
(see times t = 170 and 210): The back front of moving pedestrians in simulations of the macro- 
scopic model is sharper compared with the simulations of the stochastic model at these times. 
This can be explained by the purely hyperbolic nature of the macroscopic PDE for this range of 
the density values. 

The a = 3 regime corresponds to stronger slowdown interactions between the two groups of 
pedestrians. We observe a considerable discrepancy between the CA and PDE models at times 
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t = 140, 170 and 210. In particular, an effect similar to blocking occurs in the Monte-Carlo 
simulations of the stochastic model; this is especially evident at time t = 170. The propagation 
of two blocked groups resembles a slow diffusion over a barrier, whereas the simulations of the 
PDE model exhibit a faster propagating wave. Therefore, the PDE model with a = 3 seems to be 
outside of the range of validity of the closure approximations used in its derivation. 

We note that the results obtained using the viscous PDE model (14. lip are quite similar to the 
inviscid ones: The diffusion simply smears the fronts but no substantially different phenomena 
have been observed. 



5.2 Fully Mixed Initial Conditions 

In the second set of numerical experiments, we simulate pedestrian movement in a periodic domain, 
which is divided into 30 sectors with 15 cells in each sector (totally N = 450 cells). We initialize 
the CA simulations with the total number of 70 pedestrians (with 35 pedestrians moving in each 
direction). The initial numbers of right- and left-moving pedestrians in each sector are specified, 
while their distribution inside each sector is purely random (uniform). In dimensional units, each 
cell is 0.466m long, each sector is 7m long, and the domain is 210m long. Other parameters in 
the CA simulations are chosen to be At = 0.005 and MC = 3000. 

The corresponding initial densities for the PDE models (14.71) and (14. lip are 

p ± (x,0) = ^ioi L<x<—L, i = l,..., 30, (5.5) 

where nf and n~ are numbers of right- and left-moving pedestrians in the ith sector, respectively. 
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The size of the computational domain is L = 210 and Ax 
PDE ( 14. lip , we take the diffusion coefficient e = 0.5. 
As in Section 15.11 we take the velocities 



1. In the simulations of the viscous 



c = lm/s, Cl = c 2 = c /a, c 3 = c /(2a), 



(5.6) 



and perform two sets of numerical simulations with velocities with a = 2 and a = 3. 

The obtained results are plotted in Figures I5.3H5.101 In the a = 2 case, the CA and PDE 
results are in a very good agreement, especially when the diffusion terms are included into the 
PDE models (Figures l5.3H5.6p . By comparing Figures 



and I5.5[ one can see that the invscid 
PDE model preserves initial pedestrian clusters longer than the CA or viscous PDE ones. The 
major qualitative difference between the CA and viscous PDE models is a "blocking" phenomenon 
observed in the stochastic simulations. When the slowdown interaction is stronger (a = 3), the 
difference between the CA and PDE simulations are more pronounced: The stochastic "blocking" 
is more severe, while the PDE models develop both "blocking" and "stop-and-go" waves (see 
Figures E/FHEHQ]). 



1 





Figure 5.3: Dynamics of right-moving pedestrians in both the CA and PDE models with the fully mixed 
initial conditions (15. 5 p and velocities (15. 6p with a = 2. 



5.3 Nonhyperbolic Regime 



In both simulations described in Sections 15.11 and 15.21 the magnitude of pedestrian densities p + 
and p~ remain smaller than the critical values for the nonhyperbolic behavior (recall that the 
hyperbolicity criterion is given by (14. 9p . 

To illustrate the nonhyperbolic behavior we choose the velocities to be (15. 6p with a = 2 and 
the initial density in the simulations of the coarse-grained PDE models to be 



p + (x,0) 



0.6, 140 < k < 210, 



0. 



otherwise, 



p-{x,0) 



0.6, 186.6 < x < 233.3, 
0, otherwise, 



(5.7) 
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Figure 5.4: Time snapshots of the solutions from Figure |5~31 The solid blue line represents the CA 
solution, while the dash-dotted red line represents the inviscid (left) and viscous (right) PDE solutions. 



that is, the initial data are in the nonhyperbolic regime. The computational domain is [0,420] 
and the simulations are conducted on a Ax = 420/1280 grid. 

The initial pedestrian distribution for the CA computations is sampled from the initial density 
( I5.7jl . The number of cells in the stochastic simulations is N = 900, the cell size is 420/900 ~ 
0.4667m, the time step is At = 0.005 and MC = 3000. 

Here, we depict only p + since p~ exhibits a similar behavior. Comparison between the stochas- 
tic and inviscid PDE simulations is shown in Figure 15.111 As one can see, the PDE solution 
develops spurious large magnitude oscillations, which demonstrates that the inviscid system (14. 7p 
is ill-posed. Figure 15.121 shows p + computed using the viscous PDE model ( 14. lip . The results 
suggest that the nonlinear diffusion present in (14. lip stabilizes the PDE solution: The bigger value 
of e = 1.5 fully supresses spurious oscillations and leads to a very good agreement between the 
CA and PDE simulations. 



6 Conclusion 

In this paper, we have developed a PDE formulation for the pedestrian dynamics in a narrow street 
or corridor. The resulting 1-D system of PDEs has been systematically derived from a microscopic 
stochastic CA model. The slowdown rules in the CA model need to be specified explicitly by 
prescribing different velocities when pedestrians moving in the opposite directions are present. 
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Figure 5.5: The same as in Figure [531 but for left-moving pedestrians. 
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Figure 5.6: Time snapshots of the solution from Figure I5.5I The solid blue line represents the CA 
solution, while the dash-dotted red line represents the inviscid (left) and viscous (right) PDE solutions. 
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Figure 5.7: Dynamics of right-moving pedestrians in both the CA and PDE models with the fully mixed 
initial conditions (15. 5 D and velocities (15. 6D with a = 3. 
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Figure 5.8: Time snapshots of the solutions from Figure l57l The solid blue line represents the CA 
solution, while the dash-dotted red line represents the inviscid (left) and viscous (right) PDE solutions. 
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Figure 5.9: The same as in Figure [57] but for left-moving pedestrians. 
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Figure 5.10: Time snapshots of the solution from Figure 15751 The solid blue line represents the CA 
solution, while the dash-dotted red line represents the inviscid (left) and viscous (right) PDE solutions. 



The condition of approximate independence of neighboring cells is essential for the derivation of 
the coarse-grained PDE model and "propagate" into the macroscopic PDE formulation resulting 
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Figure 5.11: Densities of the right-moving pedestrians computed by the CA (left) and inviscid PDE 
(right) models with the initial conditions (15.70 and velocities ( 15.60 with a = 2. 



in a particular form of the flux. In addition, we have also demonstrated that the resulting PDE is 
only conditionally hyperbolic. To stabilize the resulting PDE system we have rigorously derived 
the nonlinear diffusion correction using the intermediate mesoscopic model. 

We have performed a number of numerical experiments comparing and contrasting the statis- 
tical behavior of the stochastic solution and its deterministic PDE counterparts. Our numerical 
experiments indicate that in the case of weak slowdown interactions, the average behavior of the 
stochastic CA model is reproduced quite well by the coarse-grained PDE model. On the other 
hand, the approximate independence assumption does not seem to hold in the case of the stronger 
slowdown regime, which results in a considerable discrepancy between the stochastic and deter- 
ministic results. In addition, our numerical experiments also indicate that the nonlinear diffusive 
corrections may lead to a significant improvement in reproducing the CA results using the PDE 
model. We would also like to point out that the PDE simulations are an order of magnitude faster 
than the corresponding CA ones. Therefore, the coarse-grained PDE systems have considerable 
numerical advantages in practical applications such as real-time predictions of the pedestrian traf- 
fic via filtering, investigation of pedestrian traffic flow on networks, etc. These and other issues 
will be investigated in subsequent papers. 
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Figure 5.12: Densities of the right-moving pedestrians computed by the CA and viscous PDE models 
with the initial conditions (15.71) and velocities (15.61) with a = 2. The solid blue line represents the 
CA solution, while the dash-dotted red line represents viscous PDE solutions with e = 0.5 (left) and 
e = 1.5 (right). 
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A Central-Upwind Scheme for Systems of Conservation 
Laws 

In this section, we briefly describe a semi-discrete second-order central-upwind scheme used for 
numerical solution of the 1-D systems ( 14.71) and (14.111) . which can be written as 



p t + F(p) x = (Q(p)p x ) x , 



(a.i; 



with p = (p + , p ) T , the flux F(p) = (f(p + )g(p ),f(p )g(p + )) T , an d the viscosity matrix 



Q{p) 



ec 



diag ((1 - p~) 2 + 2a ip -(l - p-) + a 3 (p~) 2 , (1 - p + f + 2a lP +(l - p+) + a 3 (p + ) 2 ) . 
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For simplicity, we divide the computational domain into finite- volume cells Cj = [Xj_i,Xj + i] of a 
uniform size Ax with x~ = j Ax and denote by 



p(x, t) dx, 



the computed cell averages of the solution p, which we assume to be available at time t. The 
semi-discrete central- upwind scheme for (1A.1|) is given by the following system of ODEs for ~pj\ 

dp j (t)_ Hj_x(t)-H J+ x(t) Pj_x(t)-P J+ x(t) 

~ir - Ax + Ax ' {A - 2) 

where the numerical hyperbolic Hj + i(t) and parabolic Pj + i(t) fluxes are constructed as follows 
(for convenience, we will omit the dependence of the computed quantities on t in the rest of the 
text). 

We first reconstruct a piecewise linear approximation, 

p(x) = pj + {p x )j{x -xj), x E Cj. (A. 3) 

The numerical derivatives (p x )j should be computed componentwise and are (at least) first-order 
approximations of p x (xj, t), calculated using a nonlinear limiter needed to avoid oscillations in the 
reconstruction (1A.3j) . In our numerical experiments, we have used the generalized minmod limiter 
(see, e.g., [H1HS]): 

where the minmod function is defined as 

{ mm j{zj}, if Zj > Vj, 
max. j{zj}, if Zj < Vj, (A.5) 
0, otherwise, 

and the parameter 9 can be used to control the amount of numerical viscosity present in the 
resulting scheme (larger values of 9 correspond to less dissipative but, in general, more oscillatory 
reconstructions) . 

The reconstruction (1A.3|) is, in general, discontinuous at the cell interfaces, x = Xj + i, where it 
has two values, which we denote by 

n ] i ■ T>, • ((>,).,■ p7-=Pj-5f(P'h (A.(i) 

These discontinuities propagate in time with one-sided local speeds that can be estimated from 
the eigenvalues A of the Jacobian matrix (14. 8[) and are calculated in the following way. We denote 
by 

R = f'(P + )g(P~) - f'(p-)g(P + ), D = [f'(p~)g( P + ) + f'(p + )g( P -)} 2 - ±f{p-)f{p + )g\p~)g\ P + ) 

and consider two possible cases: 
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If both Dj > and Dj^_ 1 > (hyperbolic regime), then 

U = \ max H + V s ?' R ^ + °} 



a . 



_I_ 

J+2 



If either < or D^\_ 1 < (nonhyperbolic regime), then 



a 



J+ i = ^ x { y/{Rj)*-Df, y/(R7 +1 )*-DY +1 } , 



The numerical fluxes are then given by 



a + 1 F(pf) - or. xFfpW.) fl + ,a: , 
= a + -a" + a+ -a" " Pj J 

P j+§ " G(p i+4 ) Ax , 2 ' 



(A.7) 



Finally, the resulting semi-discretization H A . 2 j) — f TAT7|) is a time-dependent ODE system, which 
should be numerically integrated using a stable ODE solver of an appropriate order. In our nu- 
merical experiments we have used the third-order strong stability preserving Runge-Kutta method 
(see [ion. 



Remark A.l Notice that the choice of one-sided local speeds in the nonhyperbolic regime is 
ad-hoc. However, it is important to point out that we have not tried to stabilize the inviscid 
PDE solution by increasing the amount of numerical viscosity: The solution has been stabilized 
by adding nonlinear diffusion terms rigorously derived from the mesoscopic formulation. 
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